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SUMMARY 

The principles of transfer function analysis have been applied to a 
passive optical heterodyne receiver to obtain the modulation transfer function 
(MTF) . MTF calculations have been performed based of an optical platform 
which is imaging vertically varying profiles at worst case shuttle orbit 
altitudes. An analysis of the derogatory effects of sampling (aliasing) and 
central obscurations on both resolution and heterodyne efficiency is given. 


INTRODUCTION 

One measure of performance of an optical imaging system is its ability to 
reproduce an object distribution with sufficient signal-to-noise ratio and 
resolution so as to make the information contained within the image useful. 
Generally, such a system may be characterized by its optical transfer function 
(OTF) or, in certain cases, by the modulation transfer function (MTF) (ref. 1) . 

For conventional imaging systems using either coherent or incoherent 
illumination, one usually assumes linearity in the Imaging process so that the 
cascading property of transfer function analysis applies (ref. 2) . Under this 
assumption, the MTF ' s of the individual subsystems (i.e., optics, detector, 
electronics, etc.) can be multiplied to give the overall transfer function. 

In this paper, the principles of transfer function analysis have been 
applied to a passive optical heterodyne receiver which is assumed to be 
imaging vertically varying spatial profiles at worst-case shuttle orbit 
altitudes. Results of the analysis show some interesting departures from the 
properties described above; namely, that the cascading property must be 
carefully applied and that optical receivers having obscurations, such as a 
Cassegrains, are not optimum for heterodyne- type detection. 


THEORETICAL ANALYSIS 
Imaging Considerations 

Consider an optical receiver which is imaging an object amplitude 
distribution as shown in figure 1. Using scalar diffraction theory, the signal 
amplitude, E g , in the detector plane, jr, is given by (ref. 3) 
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where h(rj is the impulse response of the imaging optics and Eg(ir) is the 
amplitude of the geometrical image of the object. The shift invariance of 
h(r) can be justified for the heterodyne applications discussed here by a 
careful examination of the various phase factors appearing in the impulse 
response function. 
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where E (r) is the local oscillator amplitude distribution in the detector 
plane, ri is the quantum efficiency, e is the electronic charge, and hv is 
the photon energy. A simple-minded classical approach is taken to obtain the 
correct expression from which the spatial frequency analysis may begin. We 
recognize that the geometrical image field, E (s,v), is a stochastic process 
which we synthesize by discrete frequency components with random phases. Now 
for a deterministic L.O. field and a quasi -monochromatic optical signal, 
equation (1) and the generalization of equation (2) combine to give 
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where < > represents an average over the ensemble of signal fields. It is 
assumed that the source, i.e., the sun, of the image field on the detector is 
spatially incoherent. The appropriate substitutions are 


<E (j3, v)E*(_s' , v) > -> Af A P (js, v ) 6 ( s_ - s_' ) 
3 3 *3 
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where P represents the image spectral radiance at the detector plane in 

2^2 2 
W/m /str/Hz and I het (f) is the current spectral power density in A /Hz. 
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Equation (3) becomes 
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where f >_ 0 and we note that the L.O. mixes with the signal field components 
at v + f and v - f. Here, we have expressed the detector overlap inte- 
gral of the L.O. field and the impulse response function as 


+oo 

1 T, ,(r)E*(r)h(r - s)d 2 r = T*(s) & h(-s) 
\ det — o — — — o — — 

—CO 


where (r) is the aperture function of the detector geometry and the 

product ? (r)E*(r) = T* (r) is simply that portion of the L.O. that is 
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transmitted by the detector aperture. 

Referring to the detector scheme of figure 2, the output current from 
the synchronous detector is 
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where H, is the total heterodyne transfer function defined by 
het • 
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and T is the optical transmission factor. The various contributions to 

H, are (1) the signal/L.O. mixing transfer function, H , representing 
het m 

carrier diffusion and transit time effects in the detector; (2) the photo- 
detector transfer function, H^, comprised of contributions due to 
capacitance, resistance and inductance; and (3) the I.F. amplifier and filter 

transfer function, H. . The square-law detector is assumed to have a unity 

ifa 

transfer function (H = 1) . The shot noise transfer function, H, , is due 

sq tr' 

only to transit time effects as opposed to H . 
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Heterodyne Transfer Function 


In this paper, we are interested in the spatial frequency response of the 
heterodyne receiver to a vertically varying object profile as shown in 
figure 1. This is different from the I.F. considerations discussed previously 
other than a knowledge of the total I.F. power. To obtain the spatial 
frequency response and, ultimately, the system modulation transfer function 
(MTF) , we assume the object scene radiance (and, consequently, the image 
scene) is linearly translated due to motion of the optical receiver, e.g. , an 
orbiting platform. This induces a translation of the image coordinates by 
an amount 
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Since the impulse response, h, is invariant, we have the output current from 
the synchronous detector as 
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Equation (4) is of the form of a convolution 
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Decomposition of I into its spatial frequency components is 
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Using the convolution theorem, we have 
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where K is the spatial frequency vector variable defined by its rectangular 

components (K ,K ) , G (K) is the object, or more specifically, the geometrical 
x y g 

image spectrum, G q (K) is the detector pupil function modulated L.O. spectrum, 
and H (K) is the coherent transfer function of the system (ref. 3). 

Equation (5) illustrates the departure of the transfer function obtained 
in a heterodyne system with that obtained in conventional imaging systems. 
Remembering that the coherent transfer function, H (K) , is equal to the pupil 
function of the optical receiver (with a suitable change in variables) (ref. 

3) , the conventional optical transfer function is proportional to 
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where G, (K) is the Fourier transform of the detector aperture function, 
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T^t (K) . In equation (5), however, we see that H (K) is modified by the 
spectrum of the L.O. /detector combination, G q (K) • The normalized convolution 

of the product G Q (K)H*(K) with its negative argument complex conjugate is 

defined as the heterodyne transfer function, G^. Functionally, then, we define 
a normalized heterodyne transfer function by 
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Heterodyne Efficiency Factor 

The denominator of equation (6) indicates that the product G o (K)H(K) 

represents the optics /L.O. detector amplitude spectrum that is transferred to 
the detector. Using Parseval's theorem, the integral 

^d 2 K| G q (K) | 2 ] H(K) | 2 is thus the power available for heterodyning out of a 
total L.O. -detector power of P q = Jd 2 ic| G q (K) | 2 . For a uniform extended 
source, we may thus define an efficiency factor 
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With this definition, the current spectrum for the synchronous detector 
(equation (5) ) becomes 


sync 


2Te 2 nV 

(K) = ■ n ~ P X G U (K) G (K) 

(hv ) 2 ° H - g - 


( 8 ) 


Equation (8) may be related to a more conventional form of heterodyne 
efficiency found in the literature (ref. 5) . The synchronous detector current 
is the inverse Fourier transform of equation (8) , i.e., 
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For a stationary scene, i.e., before translation of the image coordinates, we 
have r = 0 so that 
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Substitution of equation (8) into equation (9) and assuming a blackbody source 
of geometrical shape factor, A^(_r) , and radiance P (r) » we have 
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where mixing occurs over cin effective bandwidth 2B centered at the L.O. 

frequency, a polarization loss factor of 0.5 is included, and A (K) is the 
Fourier transform of A g(£.) • The integral portion of equation (So) has the 
form of a throughput, i.e., that portion of the image passed by the heterodyne 


500 



transfer function. The product of x and this integral is an efficiency 
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If now, we define a shot noise level due to the L.O. by 
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as the heterodyne quantum efficiency and b S j F i s t ^ le effective shot noise 

bandwidth. Note that for a uniform extended source, (K) = 6 (K) and 

equation (11) reduces to X^ et = X and ^het = OX- l n this case, the 

efficiency factor, x» which we have defined in equation (7) is equivalent 

(to within the D.C. quantum efficiency, n) to the heterodyne quantum 

efficiency, n, found in the literature (ref. 5). 
het 


System Transfer Function 

Results from the previous section may now be used to calculate the 
system transfer function, including the low-pass filter (see figure 2) , for 
the specific case of imaging a one-dimensional object through an optical 
receiver which has rectangular symmetry. This case has some physical 
significance since the resolution elements of interest in an orbiting 
heterodyne receiver are vertically varying stratospheric layers. In addition, 
to avoid scaling difficulties in the calculations we will use angular 
coordinates defined by (see figure 1) 
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where 9_ is the geometrical instantaneous field of view (I.F.O.V.) of the 
. F . 

optical receiver. 

Using equations (12) and the one-dimensional geometry, equation (6) 
becomes 
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Equation (13) assumes a plane wave local oscillator incident of the 
detector so that the detector/L.O. transfer function becomes simply the 
Fourier transform of the detector aperture. Further, the coherent transfer 
function for the optics is the pupil function (rectangular in shape) having 
a coherent cut-off frequency of K c = D^/2A, where is the diameter of the 

receiver aperture and A the wavelength. This convolution process is shown 
in figure 3. 


Equation (13) along with equation (8) gives the system transfer function 
up to the low-pass filter. Expressing G (K ;K , 0 ) and x( K 
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Finally, inclusion of the low-pass filter transfer function, H , (i.e., 

XjP 

the integration time, t) which, in this case is modeled as a running mean 
integrator, we have 


h lp ( V = sinc( ir V 

where v is the vertical component of the orbital velocity, z the receiver- 
object distance, and t the integration time. The total transfer function is 
then the modulation transfer function 
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Equations (14) and (15) can be evaluated in terms of tabulated functions 
yielding the following relations which will be used for computational purposes 
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RESULTS AND DISCUSSION 
MTF Calculations 

Equations (16) and (17) may now be evaluated for some specific parameter 
values which are applicable to the optical receiver in a space-lab type of 
scenario. A worst-case set of orbital values would be for the receiver platform 
to be at an orbital height of R = 400 Km and a tangent height, H T , of 10 Km. 

At these values, we assume that the receiver is operating in a solar occulta- 
tion mode where the sunrise or sunset velocity due to orbital motion is 
v = 2 Km/sec. The MTF and x calculations (equations (16) and (17) ) will be 
done for an I.F.O.V. of 0 = 0.5 x 10~3 rad , an equivalent optical receiver 

r 

aperture of D = 2.0", and values of integration time of t = 0.2 sec , and 

0.4 sec. Further, the value of D^ = 2.0" at a wavelength of X = 11.152 pm 

(HNO^ line) corresponds to an optics cut-off frequency of K c = D^/2X = 2278 

cycles/radian. These parameters are compatible with the values for an LHS type 
experiment using a tunable diode laser as the L.O. and associated optics for 
coupling this type radiation to a detector having the required time-frequency 
response (ref. 7) . 

The calculations are shown in figure 4. It can be shown that, for values 
of x greater than roughly 0.2 sec., the optical (heterodyne) transfer function 
dominates the MTF; and, for x somewhat less than 0.4 sec., the low pass 
filter is the dominant frequency limiting factor. Note that the angular 
frequency values can be converted to linear spatial frequency (cycles/Km) by the 
relationships of equation (12) by appropriately scaling image and object space 
by the ratio of image distance, d^ , to object distance, z. For the orbital 
values assumed, z = 2262 Km and consequently a value of 2262 cycles/rad 
corresponds to an object spatial frequency of 1 cycle/Km. Examination of the 
MTF curves shows that resolutions of the order of 1. 5-2.0 Km may be expected 
for the various integration times. 


Efficiency Calculations 

The efficiency factor (heterodyne efficiency) given by equation (17) is 

shown plotted in figure 5 for the case of the plane-wave L.O. . Two geometries 

are shown: rectangular optics (as has been previously assumed) and circular 

optics adjusted for equal optics and detector areas. The significance of the 

abscissa (2it0 K ) relative to heterodyne efficiency becomes apparent when it 
t c 

is noted that at the value of 2 tt9 f K c = 7.7 corresponds approximately to an 

image (sun) size filling the detector of one Airy Disk of the receiver 
aperture. In this region, the efficiency is in excess of 80%. 
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Sampling Error 

The calculations shown plotted in figure 5 do not include any sampling 
errors which may occur. Suppose we sample the output of the low-pass filter, 
which has been modeled as a running mean integrator, at a rate of the inverse 
of the integration time. This is equivalent to a sampled mean integration 
scheme. Under this constraint, it may be shown that for certain values of x 
the signal is undersampled. This results in an aliasing or foldover error 
which can be significant relative to the desired signal. For example, shown in 
figure 6 is the total MTF for the values of orbital and system parameters 
previously stated. Two integration times are considered: x = 0.2 sec and 

x = 0.4 sec. If we define the sampling error as the ratio of the "foldover" 
amplitude on the MTF plot to the amplitude of the MTF itself, i.e., a white 
signal spectrum, we see that the error for 0.4 sec. is approximately 40% at 
0.5 cycle/Km frequency and considerably worse for higher values of spatial 
frequency. Conversely, for x = 0.2 sec and the correspondingly higher 
sampling rate, the sampling error is negligible. 


Heterodyne Receivers With Obscurations 


Telescopes having central obscurations such as Cassegrains are often used 
for imaging a source. If this type receiver is used as a collector for 
heterodyne-type detection, one needs to compare the efficiency, x» and the 
heterodyne transfer function, , with that obtained for the unobscured case. 


In figure 7, we consider the effects of receiver apertures having obscura- 
tion ratios of 0 and 20% for 0 = 0.2 x 10~ 3 and 0.5 x 10 -3 rad. Note the 

£ 

enhancement of response in the 2000 cycles/rad region at the expense of that 

near 1000 cycles/rad for 20% obscuration and 0 = 0.5 x 10 -3 rad. The 

£ 

effects of obscurations are more pronounced for square as opposed to circular 
geometries. An unobscured conventional MTF discussed earlier is plotted for 
0 p = 0.5 x lo~3 rac j showing a somewhat reduced frequency response 
characteristic from the heterodyne MTF. For a smaller detector 
(0,-, = 0.2 x lo -3 rad ), a 20% obscuration tends to assume the shape of a 

r 

conventional MTF. 


In figure 5, we assume a receiver aperture having obscuration ratios of 

20% and 50%. For the values of 0_ and K used earlier, X = 2u0 K =7.2, 

F c F c ' 

and comparison of the various efficiency curves at this value shows striking 

differences. For the 50% case, one sees that the heterodyne efficiency is 

virtually zero while for the 20% case a relative efficiency of slightly 

greater than 0.2 is achieved. This compares with a value of greater than 0.8 

in the unobscured case. Note further there is a "peaking" of the efficiency 

curves for obscured systems. The rule-of-thumb requirement of one Airy disk 

on the detector for "good" efficiency no longer holds but rather the source 

image needs to be less than this value to achieve the maximum efficiency for 

that particular system. The effect can be explained on the basis of the 

overlap integral (equation (2) ) of the L.O. field and signal field 
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distributions. Thus, the 'diffracted field due to the central obscuration is 
out of phase with that of the primary diffracted field and, as the size of the 
detector and/or optics increases, the cancellation tends to be more complete. 


CONCLUDING REMARKS 

The analysis of a passive heterodyne receiver with respect to its imaging 
performance (transfer function) and its heterodyne efficiency shows some 
interesting departures from the results which are obtained in strictly coherent 
or incoherent imaging systems. For example, the cascading property of MTF 
analysis must be carefully applied since the coherent transfer function of the 
optical receiver and that due to the L.O. -detector combination are not 
separable but are related by the convolution of their products. Application 
of these results to the specific case of a space-lab type optical heterodyne 
receiver (LHS) shows that resolutions of the order of 1. 5-2.0 Km are possible 
for worst-case type orbital scenarios. 

Further, comparison of obscured-type receivers (e.g., Cassegrains) with 
unobscured receivers shows that both resolution and efficiency are severely 
degraded in an obscured-type receiver and consequently should not be used for 
a passive heterodyne detection scheme. 
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Figure 2.- System transfer functions. 
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Figure 5.- Efficiency factor, X/ versus system parameter 2TT0.pK c for various 

receiver geometries. 



Figure 6.- Aliasing error for 2.5 Hz sampling rate at worst-case shuttle 

orbit. 



Figure 7.- Comparison of transfer functions for various geometries. 

curve represents conventional MTF with 0.5 mrad I.F.O.V. 
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